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We perform measurements of the third moment of atom number fluctuations in small sHces of a 
very elongated weakly interacting degenerate Bose gas. We find a positive skewness of the atom 
number distribution in the ideal gas regime and a reduced skewness compatible with zero in the 
quasi-condensate regime. For our parameters, the third moment is a thermodynamic quantity whose 
measurement constitutes a sensitive test of the equation of state and our results are in agreement 
with a modified Yang- Yang thermodynamic prediction. Moreover, we show that the measured 
skewness reveals the presence of true three body correlations in the system. 
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Measurements of higher-order correlations and the 
density fluctuations, in particular, are becoming an in- 
creasingly important tool in the studies of ultracold quan- 
tum gases. Such measurements are able to probe quan- 
tum many-body states of interacting systems, often giv- 
ing access to key quantities that characterize the sys- 
tem. This is particularly true for one-dimensional (ID) 
gases, where the effects of fluctuations are enhanced com- 
pared to 3D systems and govern the rich underlying 
physics. Zero-distance (or local) second- and third-order 
correlation functions have been probed in several ultra- 
cold gas experiments by measuring the rates of two- and 
three-body inelastic processes such as photoassociation 
and three-body recombination[lH3|. Such measurements 
enabled the study of the strongly correlated regime of 
'fermionization' in a ID Bose gas. 

An alternative experimental technique is the in situ 
measurement of atom number fluctuations in a small de- 
tection volume, achievable through the analysis of noise 
in absorption images. The fluctuation variance (or sec- 
ond moment) provides information about an integrated 
non-local density-density correlation function. In addi- 
tion, under adequate experimental conditions, the vari- 
ance can render as a thermodynamic quantity and there- 
fore such measurements constitute a probe of the ther- 
modynamic properties of quantum gases [a; fi L alterna- 
tive to the analysis of density profiles [2, LD, l9| or mo- 
mentum distributions in the context of Tan's relations 
[lOJ. In situ measurements of atom number fluctua- 
tions in weakly interacting quasi- ID Bose gases were 
used to probe the crossover from the nearly ideal gas 
regime, where bosonic bunching is present, to the quasi- 
condensate regime, where the density fluctuations are 
suppressed [J, |6|. In fermionic systems, sub shot-noise 
atom number fluctuations were observed in a degenerate 
Fermi gas [ll|,ll2|- Combining this with the measurement 
of compressibility of the gas deduced from the known den- 
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FIG. 1. (Color online) (a) Scheme of the imaging setup. The 
probe laser crosses the atomic cloud (red dot) before its reflec- 
tion from the chip surface and detection on a CCD camera, 
(b) Typical m-situ absorption image; the scale on the colorbar 
shows the optical density. The pixel size in the object plane 
is A = 4.5 /xm. (c) Typical longitudinal density profile (solid 
curve), together with the mean profile (dashed red curve). 



sity profile and confining potential, such measurements 
have been shown to provide reliable thermometry |ll| . 

In this paper we expand the arsenal of probes of higher- 
order correlations in quantum gases by measuring the 
third moment of atom number fluctuations. This is done 
using in situ absorption imaging of an ultracold gas on 
an atom-chip setup sketched in Fig. [T] (a). We probe a 
weakly interacting quasi-lD Bose gas. We have measured 
a positive third moment of the atom number distribution 
in a degenerate gas within the ideal gas regime and within 
the crossover towards a quasi-condensate. In the quasi- 
condensate regime the measured third moment is com- 
patible with zero. The third moment of the atom number 
distribution is linked to the third-order correlation func- 
tion and our measurements demonstrate the presence of 
true three-body correlations in the gas. Apart from this, 
we show that the measured third moment is related to 
a thermodynamic relation that involves a second-order 
derivative, and therefore the technique can be used as a 
sensitive probe of the thermodynamics of a quantum gas. 



Our quasi-lD Bose gases are produced using ^^Rb 
atoms in the hyperfine state \F = 2,m = 2). A very 
elongated loffe magnetic trap with a longitudinal oscilla- 
tion frequency ranging from 5.0 Hz to 8 Hz and a trans- 
verse oscillation frequency uj±/2Tr ranging from 3 kHz to 
4 kHz is realized using on-chip micro-wires and an exter- 
nal homogeneous magnetic field. Using forced rf evapora- 
tion, we produce ultracold clouds at temperatures from 
T = 20 nK to 500 nK. The longitudinal rms size L of 
the cloud ranges from ~50 ^m to ~100 ^m. As shown 
previously [J|, under these conditions such gases explore 
the crossover from the ideal gas regime to the quasi- 
condensate regime, and the underlying physics lies in the 
ID regime or in the crossover from ID to 3D [6|. 

In situ measurements of density fluctuations are per- 
formed using absorption images such as the one shown 
in Fig. [T] (b). The details of our imaging and calibra- 
tion techniques are described in EPAPS Ref. [13|. As the 
transverse size of the trapped cloud (< 500 nm rms) is 
much smaller than the pixel size (4.5 /im), the only in- 
formation in the transverse direction is the diffractional 
and motional blur on the image. By summing the atom 
number over transverse pixels, we reduce the notion of a 
pixel to a segment of length A and derive from each image 
the longitudinal density profile [Fig. [1] (c)]. We perform 
a statistical analysis of hundreds of images taken under 
the same experimental conditions [J, |6| . For each pro- 
file and pixel we extract 5N = iV — (iV), where (N) is 
given by the average density profile. To remove the ef- 
fect of shot-to-shot variations in the total atom number 
A^tot, the profiles are ordered according to A^tot and we 
use a running average of about 20 profiles. As will be 
explained below, the longitudinal confining potential is 
irrelevant and each SN is binned according to the corre- 
sponding mean atom number in the pixel (N). For each 
bin, we compute the second and third moment of atom 
number fluctuations, {SN"^) and {5N^). The contribution 
of the optical shot noise to these quantities is subtracted, 
although it is negligible for (SN^). 

The measured third moment of the atom number fluc- 
tuations, {SN^)m, is plotted in Fig. [2] for two different 
temperatures. For the higher temperature [Fig. [2] (a)], 
we observe a positive value of {SN^)m that increases with 
the average atom number (A^). At a smaller temperature 
[Fig.[2](b)], {5N^)m initially grows with (A^) and reaches 
a maximum, before taking a value compatible with zero 
at large ( A^) . The corresponding second moments or vari- 
ances {SN'^)m are shown in the insets. A finite third mo- 
ment indicates an asymmetry of the atom number distri- 
bution, which is usually quantified by the skewness of the 
distribution, s,„ = {5N^)m/ {5N'^)m , shown in Figs. [2] 
(c) and (d). Before discussing the physics and theoreti- 
cal understanding of these results, we first describe how 
the measured moments {5N^)m and {5N'^)m are related 
to the true moments (SN^) and {SN'^). 

The measurements of atom number fiuctuations are 
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FIG. 2. Measured third moment (open circles) of the atom 
number fluctuations versus the mean atom number per pixel, 
for temperatures of 376 nK (a) and 96 nK (b). The insets 
show the corresponding atom number variances. The error- 
bars are the statistical errors. Graphs (c) and (d) show the 
skewness Sm corresponding to (a) and (b), respectively. The 
theoretical predictions, scaled by K2 = 0.55 and K3 = 0.34 for 
(a) and (c), and by ^2 = 0.52 and ^3 — 0.31 for (b) and (d), 
are shown for comparison: solid lines - the modified Yang- 
Yang prediction; dashed lines - the ideal Bose gas prediction; 
dash-dotted lines on (b) and (d) - the quasi-condensate pre- 
diction; dotted lines - the shot noise limit (A^) . 



affected by the finite spatial resolution due to both the 
optical resolution and the diffusion of atoms during the 
optical pulse, which cause the absorption signal from 
each atom to spread over several pixels and blur the 
image. Denoting by A the impulse response function 
of the imaging system, the impulse response for the 
pixel [0, A] is J-{zq) ~ L dzA{z — zq), and the mea- 
sured atom number fluctuation in the pixel is given by 
6Nm — /_ dzQj-{zQ)5n{zQ), where (Sn(zo) is the local 
density fluctuation. For the parameters explored in this 
paper, the expected correlation length l^ of density fluc- 
tuations [IJI is smaller than 0.5 /im. This is sufficiently 
smaller than the width of A so that we can assume that 
the density fiuctuations have zero range. Moreover, since 
the resolution and the pixel size are much smaller than 
the longitudinal size of the atomic cloud, we can assume 
that the gas is locally homogeneous with respect to z. 
Then, the measured second and third moments of den- 
sity fiuctuations can be obtained as 

{5N^)„, = {5N^) j:^^dzo Hzof/A = «2 (SN^), (1) 
{SN'),-a - {6N') /+^d2o -F(;^o) V A - «3 (SN^), (2) 

where {SN'^) and {SN^) are the respective true moments, 
whereas K2 and K3 are the reduction factors. For low 
enough linear densities, the gas lies in the nondegenerate 
ideal gas regime. Then the fluctuations are almost that of 
a Poissonian distribution, so that (SN^) ~ {5N^) ~ (A^), 
and the reduction factors may be deduced from a linear 



fit of the measured fluctuations versus (TV) , wiiere (TV) is 
experimentally determined absolutely. However, such a 
deduction is difficult in very cold clouds where only few 
pixels lie in the nondegenerate ideal gas regime. 

We thus develop an alternative method that uses the 
measurement of the atom number correlation 

^i^i+j— /!;at2\ ~ r+oo, ^, so ^'^) 
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between the pixel i and the adjacent (j = l) or the next- 
neighbor (j = 2) pixels. Such correlation arises due to the 
contribution of an atom to the absorption in both pix- 
els. Making a Gaussian ansatz for the impulse response 
function A^ the rms width of A can obtained by fitting 
Eq. ^ to the measured correlations Q^+i [13|. The re- 
duction factors K2,3 can then be deduced from Eqs. ([T]) 
and ^, resulting typically in K3 ~ 0.3 and K2 — 0.5 for 
our data. The result is in good agreement with the slope 
of {SN^rn at small (TV) [see the inset of Fig. [2] (a)]. 

Turning to the discussion of the physics behind our ex- 
perimental results, we first point out that the third mo- 
ment of atom number fluctuations is actually a thermo- 
dynamic quantity when, as in our experiment, the pixel 
size is both much larger than the characteristic corre- 
lation length of density fluctuations l^ and much smaller 
than the cloud length L, /^ ^ A <C -L. Then a local density 
approximation is valid and the gas contained in a pixel 
can be well described by a grand canonical ensemble, in 
which the rest of the cloud is acting as a reservoir that 
flxes the chemical potential fi and the temperature T. 
Denoting by Z the grand-canonical partition function, we 
have (N) = {kBT/Z)dZ/dn, (N^) ^ {klT^/Z)d'^Z/dfi'^ 
and (N^) = {k%T^/Z)d^Z/dfi^. From the first two equa- 
tions, we obtain the well known thermodynamics relation 
((5TV^) =kBTd{N)/d^, whereas the three equations give 
the following relation 
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(4) 



where (N) = nA, and n is the linear density of a gas ho- 
mogeneous along z. Thus, the knowledge of the equation 
of state (EoS) n = n{ii, T) is sufficient to predict the 
third moment of the atom number distribution. Note 
that a more traditional form of the EoS for pressure P 
can be readily deduced from n(^, T) using Gibbs-Duhem 
relation. n={dP/dii)T, leading to P = J!^ n{fi' ,T)dii' . 
We now compare our measurements with the predic- 
tions from different models for the EoS n{fi, T). The tem- 
perature of the cloud for the case of Fig. [2] (a) is deduced 
from an ideal Bose gas fit to the wings of the density 
profile [6|. For the data of Fig. [2] (b), corresponding to 
the quasi-condensate regime, such wings are vanishingly 
small and hard to detect. In this case we deduce the 
temperature [6| from the measurement of density fiuc- 
tuations in the cloud centre using the thermodynamic 



relation {SN^ =k%T'^d{N)/dfi and the EoS of a quasi- 
condensate (see below). 

The predictions from the equation of state for an ideal 
Bose gas are shown by the dashed lines in Fig. [5J For a 
highly nondegenerate (or classical) gas, corresponding to 
small (TV), this model predicts {SN^) ~ {SN"^) ~ {N) as 
expected for a gas of uncorrelated particles. When the 
gas becomes degenerate with the increase of (N) , the con- 
tribution of the quantum-statistical exchange interaction 
term to {6N^) is no longer negligible, and (SN^) becomes 
larger than the shot- noise term (N) . Such an increase is 
observed in the experimental data in Fig. [2] (a). However, 
the ideal Bose gas model strongly overestimates the third 
moment with further increase of (N) and we eventually 
observe large discrepancy between the predictions of this 
model and the experimental data. The discrepancy is due 
to the repulsive interactions between the atoms, which 
reduce the energetically costly density fluctuations. 

Describing the effects of atomic interactions beyond 
the perturbative regime is a challenging theoretical prob- 
lem. However, a purely ID Bose gas with contact repul- 
sive interactions is a particular case of a many-body prob- 
lem for which an exact solution is available through the 
Yang- Yang thermodynamic formalism [l5| in the entire 
parameter space. For the temperatures corresponding to 
Fig. [2] (a)-(b), the ratios of kBT/niu± are 2.6 and 0.50, 
respectively, implying that the population of the trans- 
verse excited states is not negligible. Accordingly, we 
use a modified Yang- Yang model [7|, in which the trans- 
verse ground state is described within the exact Yang- 
Yang theory [16|, whereas the transverse excited states 
are treated as ideal ID Bose gases. This model has been 
shown to be valid for our parameters until the quasi- 
condensate regime is reached [6| . The corresponding pre- 
dictions are plotted in Fig. [5] and show that the model 
accounts for the measured (SN^) very well. 

In the quasi-condensate regime [corresponding to 
(N) > 70 in Figj2] (b)!, where the density fluctuations 
are suppressed [171 ll8| i the EoS can be obtained numer- 
ically from the 3D Gross-Pitaevskii equation and is well 
described by the heuristic function fj, = huj j_{y'l + 4na — 
1) |19|. In contrast to the modified Yang- Yang model, 
this EoS accounts for the transverse swelling of the cloud 
due to interatomic interactions and better describes the 
measured variance [see the inset of Fig. [5] (b)] . The mea- 
sured third moment is compatible with this EoS. 

To unveil the role of many-body correlations, which un- 
derly the measured density fluctuations while remaining 
hidden in the thermodynamic analysis, we consider the 
ID two- and three-body (fc = 2,3) correlation functions, 

5W(zi, . . . , Zk) = {^PHzi) . . . i^^ZkMzk) . . . ^(zl))/r^^ 

(5) 
where ip{z) = jdxdy ip{x,y,z) and ip {ip^) is the bosonic 
field annihilation (creation) operator. Using standard 
commutation relations and the expression (TV^) = (iV)^ + 
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FIG. 3. The measured three-body integral Hm versus (A'^), 
corresponding to the data of Figs. [2] (a) and (b), respectively. 
The solid (dash-dotted) curves are the thermodynamic predic- 
tions from the modified Yang- Yang (quasi-condensate) model. 



{N) + n'^ JJdzidz2[g'^'^Hzi,Z2) - 1], we find 
{SN^) = {N) + 71^ JJJdz,dz2dz3 [g^^\zuZ2,Z3) - l] 
-3{N)n'JJdzidz2[g'^^\zi,z2) - 1] 
+3n'^JJdzidz2[g^^\zi,Z2) - 1], (6) 

where the integrals are in the interval [0, A]. As we see 
the third moment of atom number distribution depends 
on both the g'-^^ and g'-'^^ functions. Moreover, g'-^^ con- 
tains a contribution from g'^-* since, when one of the three 
particles is far from the other two, g^^^ reduces to g^^^. To 
remove such contributions, we introduce the /i-function 

h{zi,Z2,Z3)=2 + g^^\zi,Z2,Z3) (7) 

-[g^'Hz,,Z2)+g'^'Hz2,Z3) + ~g<^'\z,,zs)], 

which is nonzero only for zi, Z2 and Z3 being all in the 
vicinity of each other. Such a decomposition has been 
previously used in the description of weakly correlated 
plasmas [20|, with the approximation /t, ~ being used 
to truncate the BBGKY hierarchy. 

Using the ^.-function, Eq. (J6|) can be rewritten as 

(6N^) = {N)+3n^JJdzidz2[g'^^Hzi,Z2) - 1] 

+ n^ JJJ dzidz2dz3h{zi, Z2, Z3). (8) 

Here, the first two terms represent one- and two-body 
effects, with the second term being equal to 3{6N'^) — 
3(iV). The contribution of true three-body correlations 
to (SN^) comes from the three-body integral 

H^{SN'^)+2{N)-3{SN^)^n^JJJdzidz2dz3h{zi,Z2,Z3). 

(9) 
In Fig. [3l we plot the measured value of "H. More 
precisely, taking into account the reduction factors K2 
and K3, we plot H = {6N^)r,Jn3 + 2(iV) - i{5N^),n/n2- 
We observe nonzero values of "H, which is a signature 
of the presence of true three-body correlations in the 
gas: "H is positive within the ideal gas regime and in the 
crossover region towards the quasi-condensate [see Fig. [3] 
(a)], whereas it is negative in the quasi-condensate regime 



[Fig. [31(b)]. The results are in agreement with the ther- 
modynamic predictions of the modified Yang- Yang and 
the quasi-condensate models. 

In summary, we have measured the third moment of 
density fluctuations in an ultracold quantum gas. This 
quantity reveals the presence of true three-body correla- 
tions in the system. Moreover, for sufficiently large pix- 
els, such measurements constitute a very sensitive probe 
of the thermodynamics of the gas. As the third moment 
is related directly to the second-order derivative of the 
equation of state n{^,,T), the method lends itself as a 
high-precision tool for discriminating between alterna- 
tive theoretical models, and can be applied to a broad 
class of ultracold atom systems. For example, intrigu- 
ing opportunities are in the understanding of the role 
of higher-order correlations in thermalisation of isolated 
quantum systems [21| and in the study of thermodynam- 
ics of more exotic many-body systems where three-body 
effects, such as Efimov resonances |22|, may lead to differ- 
ent signatures in the second- and third-order correlations. 
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03 and by the Australian Research Council. 
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FIG. 4. (Color online) Typical absorption images used for cal- 
ibration and the analysis of density fluctuations: (a) - taken 
with a nearly resonant probe, after a small time of flight of 
0.5 ms [same as Fig. 1 (b) of the main text, except with a 
larger field of view]; (b) - taken with a probe detuned by 
5 MHz and a time of flight of 1.5 ms; (d) - taken with a 
resonant probe and a time of flight of 2.2 ms to measure the 
total atom number. The scales on the colorbars correspond 
to optical densities. 



Precise measurement of the atomic density profile 
and determination of imaging resolution 

A probe laser beam, locked onto the D2 transition at 
the wavelength A = 780 nm, is reflected from the chip 
surface (covered by a gold mirror) after passing through 
the atomic cloud. The shadow image of the atomic cloud 
is then recorded on a CCD camera, with quantum effi- 
ciency larger than 90%. The diffraction-limited optical 
resolution has an rms width of 1.0 /xm. Great care has 
to be taken in absorption imaging to allow for a reliable 
measurement of the atom number in each pixel. 

We have chosen a configuration that maximizes the ab- 
sorption efficiency. To achieve this, we focus the probe 
beam onto the chip surface using a cylindrical lens as de- 
picted in Fig. 1 (a) of the main text. The 1/e^ size of 
the beam in the focused direction x is 50 /xm, which is 
smaller than the distance of the atomic cloud from the 
chip surface so that the beam crosses the atomic cloud 
only on its way to the chip and not after its reflection. 
Images as in Fig. 3] (a)-(c) are taken after switching off 
the wire currents so that only the external homogeneous 
magnetic field, whose orientation is close to the y direc- 
tion [see Fig. 1 (a) of the main text for the axis defini- 
tion[. remains switched on. With this geometry, using 
a cr+-polarized probe beam we address only the closed 
transition |F = 2,m = 2) -> \F' = 3, m' = 3), and the 
absorption cross-section at low intensity takes its maxi- 
mum value CTo = 3A^/27r. 

The absorption is measured by taking two pictures, 
the first one with the atomic cloud present and the sec- 




60 80 100 120 
Nbl 



FIG. 5. Experimental determination of the correction to the 
Beer-Lambert law. The solid curve is a fit to the experimen- 
tal data with a third-order polynomial, the straight dashed 
line being the linear contribution. Nbl is obtained by using 
the Beer-Lambert law and summing over transverse pixels in 
images taken with a nearly resonant probe and a small time 
of fiight of 0.5 ms [see Fig. |4](a)]. A'^ is obtained from images 
taken with a probe detuned by 5 MHz and a time of flight 
of 1.5 ms, for which the Beer-Lambert law is adequate [see 

Fig.i](b)[. 



ond one without the atoms. The atom number in a 
given pixel Np is estimated from the Beer-Lambert law 
Nj, = ln(7V2/7Vi)A2/CT, where Ni and N2 are the photon 
numbers in the pixel on the first and second image, re- 
spectively, and a is the atomic absorption cross-section. 
With the transverse size of the cloud being smaller than 
the pixel size, no information is available in the transverse 
direction, and the number of atoms Nbl in an effective 
pixel of size A can be obtained by summing Np over of 
the transverse pixels. However, as the use of a resonant 
probe at high atomic densities produces high optical den- 
sities (up to 1.5) this naive procedure fails to correctly 
estimate the true atom number N in the pixel. Firstly, 
when the transverse extension of the cloud is smaller than 
both the pixel size and the optical resolution, the Beer- 
Lambert law underestimates the true atom number due 
to the concavity of the logarithm function, as already 
pointed out in [l|. Moreover, the validity of the Beer- 
Lambert law is questionable for high atomic densities due 
to nontrivial reabsorption effects that may arise [2| . In or- 
der to reduce these effects, without decreasing too much 
the absorption, we use a near resonant probe and enable 
the cloud to spread transversally during a small time of 
flight of about 0.5 ms - sufficient to reduce the effects 
of high atomic densities, but small enough so that the 
atom number fluctuations in a pixel are barely affected. 
The measured atom number Nbl, however, still deviates 
from the true atom number A'^ and we introduce a func- 
tion / deflned as iV = I{Nbl) to describe the deviation 
from linearity at high optical densities. 

The function / is deduced from the comparison, in 
each effective pixel, of the measured Nbl with the cor- 
rect atom number N . The latter itself is measured as fol- 
lows. First, the correct proflle is obtained using images 
[such as that shown in Fig. S] (b)] taken with a detuning 
of about 5 MHz that reduces the absorption cross-section 



and a time of flight of ~ 1.5 ms that permits a transverse 
expansion of the cloud. We checked experimentally that 
these parameters ensure the validity of the Beer-Lambert 
law, while the expansion is small enough as to retain the 
longitudinal profile essentially unaffected. Second, the 
absolute normalization of the atomic density profile (or, 
equivalently, a measure of the absorption cross-section 
a) , is deduced form the knowledge of the total atom num- 
ber. The latter is measured using a resonant probe with a 
time of flight of ^2.5 ms, as in Fig.[4](c). With such time 
of flight, the cloud transverse expansion is sufficiently 
large and the atomic density is small enough as to render 
the Beer-Lambert law applicable. Atomic saturation is 
taken into account via the formula a = (To(1 -|- ///sat), 
where / is the intensity of the probe beam and Igat is the 
saturation intensity. A fit of the measured absorption 
versus / gives /gat = 1-4(1) mW/cm^, which is close to 
the reported value of 1.62 mW/cm^ [3|. The remaining 
discrepancy could be because of possible underestimation 
of the intensity of the probe beam seen by the atoms due 
to the losses during the reflection of the beam from the 
gold mirror and during the transmission through the op- 
tical lenses. Finally, the function / is estimated by fitting 
the experimental points N versus Nbl with a third-order 
polynomial; as shown in Fig. [5] 

All these calibrations are performed using values aver- 
aged over tens of experimental realizations. The images 
used for the analysis of fluctuations [as in Fig. |4] (a)] and 
those used for calibration [as in Figs. |4] (b) and (c)] are 
taken in an alternated way - typically one picture of type 
(b) and then (c) after every 15 images of type (a) - to 
eliminate the dependence on noise arising from long-time 
magnetic field and thermal drifts. 

The normalization procedure described above, while 
compensating for the effect of the small transverse size 
of the atomic cloud, does not compensate for a possible 
error induced by short-scale longitudinal density fluctua- 
tions. However, in our experimental situation, those fluc- 
tuations are smeared out by the atomic diffusion during 
the probe pulse and are small. 
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FIG. 6. Correlation of atom number fluctuations between 
adjacent (open squares) and next-neighbor (crosses) pixels. 
The solid lines are the predictions for an impulse response 
function with an rms width of 2.0 /im, obtained by fitting 

The precise calibration of the atom number measure- 
ment described above is not sufficient for the analysis 
of atom number fluctuations. Indeed, the fluctuations 
are affected by a finite spatial resolution due to both 
the optical resolution and the diffusion of atoms during 
the optical pulse, which cause the absorption signal from 
each atom to spread over several pixels and blur the im- 
age (see main text). The finite spatial resolution is also 
responsible for a nonvanishing correlation between the 
atom numbers measured in nearby pixels. In fact, we 
make use of this correlation for experimental determina- 
tion of the spatial resolution of our imaging system. In 
Fig. [5] we show the experimental data for the correlation 
coefficients Ci^i+j [see Eq. (3) of the main text[, corre- 
sponding to the atom number fluctuations in the adjacent 
(j = 1) and next-neighboring (j = 2) pixels, for the ex- 
perimental data of Fig. 2 (c) of the main text. Fitting 
these correlation coefficients, we extract the rms width 
5 of the impulse response function A and we find that 
(5 = 2.0 yum for this data set. 
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